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ABSTRACT 


This study presents two methods of analysis for the 
supersonic oscillating cascade with subsonic leading edge. 
A relatively straightforward solution is developed for the 
slowly oscillating finite and infinite flat-plate cascades 
which provides simple analytical expressions for the unsteady 
pressure distributions. Comparison with other solutions 
is generally excellent. Some additional linear topics 
including resonance and a unique inflow condition are also 
treated. In addition a nonlinear method of characteristics 
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l. INTRODUCTION 


l.l General. Unsteady flow in turbomachines can be broadly 
divided into two categories: circumferential asymmetry in 


the flow through the turbomachine and flow perturbation due 
to the self-excited oscillatory motion (flutter) of the 
internal blading [1]. As shown in figure 1.1.1, taken from 
[2], at least four types of flutter may occur in a compressor 
rotor. The two stalled-flutter phenomena are encountered 
when compressor or fan blading is subjected to incidence 
angles near or exceeding the stall angle, and so are off- 
design phenomena. Subsonic stalled flutter has long plagued 
the engine designer [3]. Choke flutter and supersonic unstalled 
flutter are lows ande instabilities occurring when the 
flow over the blade is, respectively, transonic along most 

of the chord, or supersonic over the outer span portions [3]. 
Supersonic unstalled flutter (hereafter referred to simply 

as supersonic flutter), which has appeared only recently as 

a significant problem in turbomachinery, is the subject of 
this investigation. Primarily a concern in compressors and 
fans [1], it is an instability born of the recent technologi- 
cal advances aimed at increasing the thrust-to-weight ratios 
of turbomachines: increased tip speeds, composite blades, 
lighter discs and blades, and the elimination of part-span 
shrouds [2,3]. Some insight into the nature of supersonic 


Mineo meuceoainedssy a consideration of its characteristics: 
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Figure 1.1.1 Compressor map showing boundaries for four 
types of flutter. 


l. Its flutter boundary imposes a high-speed limit on 
compressor operation [2]. 

2. It is a low-incidence, attached-flow instability that 
eeaeoeccur ae the design condition (figure 1.1.1). 

5. Flutter amplitude decreases with increasing back 
pressure on the rotor [2]. (For this reason the present 
analysis considers the most critical case of zero pressure 
rise across the rotor.) 

4. . At the flutter condition all the blades are observed 
to flutter at the same frequency with a constant interblade 


phase angle [2]. 





5, Both single-degree-of-freedom torsional flutter (pitch) 
and two-degree-of-freedom torsional-flexural flutter (pitch 
and plunge) can occur [3,4]: 

6. The flutter boundaries are significantly affected by 
airfoil shape [4,5]. 

7. The primary flutter mechanism is an inviscid phenomenon. 
ME the unsteady aerodynamic force lags the blade motion 
by a phase angle which is sufficient to ensure positive 
work over the entire cycle, then flutter can occur Rd 

It is clear that reliable techniques for predicting the 
unsteady pressure distribution on a fluttering blade are 
required if one hopes to predict and avoid the onset of 
flutter. With this end in mind the cascade has been exten- 
sively employed as a theoretical abstraction of a compressor 
rotor. The unsteady flow is considered in an annulus of 
differential radial height within which equispaced airfoils 
M radially disposed. Unwrapping such an annulus gives 
rise to the two-dimensional approximation of plane flow 
through a cascade of airfoils [1]. The cascade is termed 
linear if the airfoils are flat plates of zero thickness. 
As seen in figure 1.1.2, a rotor does not require a super- 
sonic axial flow to have a supersonic velocity component 


relative to, and parallel to, the blade axis. Indeed this 


D ————— 


Е. Агпо1ді её а1. іп "Supersonic Chordwise Bending 
Flutter in Cascades," Pratt € Whitney Report No. PWA-5271 
of 31 May 1975 show that choráwise bending can also occur. 
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Figure 1.1.2 Supersonic cascade classification. 





distinction provides a convenient method of classifying 
supersonic cascades. If the axial velocity is subsonic, 
the cascade is said to have a subsonic leading edge; in the 
other case it has a supersonic leading edge [6]. The physi- 
Seemotenificance of this distinction is clear from figure 
1.1.2. In the supersonic-leading-edge cascade each blade 
influences only its immediate neighbors. When the cascade 
has a subsonic leading edge, however, disturbances generated 
by the upstream portions Of each blade affect the flow around 
each subsequent blade in the cascade; furthermore the flow- 
field at the underside of each blade is influenced by the 
wake disturbances of all preceding blades. Since the supersonic- 
leading-edge cascade has been EE investigated [7,8,9, 
10], the more' difficult subsonic-leading-edge case, which is 
also the more interesting in practice, is considered in the 
present study. Another cascade classification is frequently 
employed: if each blade is assumed to be preceded and followed 
by infinitely many others, the cascade is said to be infinite; 
if there is an identifiable first blade — a necessary pre- 
requisite for many numerical solutions — the cascade is 
termed finite. 

A considerable amount of effort has been expended in 
attempting to predict — numerically and analytically — 
the unsteady pressure distribution over a blade in an oscilla- 
ting supersonic cascade with subsonic leading edge. A 
brief discussion of some of the more noteworthy advances 


follows. 





BD The Linear SUBsontes Lcda Edge Cascade. Verdon 


and McCune [11] formulated a linear boundary value problem 
for the subsonic-leading-edge configuration, solving it 

for velocity potential. This permitted calculation of the 
pressure along the entire upper surface of an arbitrary blade 
and along the lower surface from the leading edge to the 
point of impingement of the trailing-edge Mach wave from 

the preceding blade. In a more recent paper Verdon [12] 
formulated a second boundary value problem, in which pressure 
is the basic dependent variable, which permitted determination 
of the pressure along the entire lower surface as well. 
Nagashima and Whitehead [13] and Goldstein [14] obtained 
solutions by replacing the cascade blades with doublet dis- 
tributions.  Kurosaka [6] introduced what he terms ane 
passage approach in applying карс transform techniques 

to a slowly-oscillating cascade, obtaining closed-form 
analytical expressions for the pressure distributions of 
certain configurations. He recently [15] extended his solu- 
ños to include arbitrary frequencies of oscillation. i 
assuming that a space periodicity of the flow exists from 

one blade passage to the next, Kurosaka claims that he 

avoids complications arising from the breakdown of linear 
acoustic theory in the far field [4]. Platzer and Brix 

[16] and Platzer and Bell [17] employed the method of char- 
acteristics to solve directly for the perturbation veloci- 
ties and pressures everywhere in the flowfield. Their 


computer program has been thoroughly checked against existing 








Figure 1.2.1 Limits of applicability of Chadwick's solution. 


solutions [11,12,13] and is extensively employed in the 
present study as a method of comparison for both the analy- 
tical and numerical solutions presented. Chadwick [5] used 
Sauer's solution [18] of the unsteady potential equation, 
accurate to a first order in frequency of oscillation, to 
determine the pressure distribution on a given blade from 
the leading edge to the point of the first reflection (point 
A in figure 1.2.1) of disturbances from the preceding blade. 
He obtained exact analytical agreement with Kurosaka [6] over 
that portion of the blade for the single-degree-of-freedom 
e сазсаде. Chadwick's technique has the decided 


advantage of being much simpler than that of Kurosaka, easily 





yielding closed-form expressions for any cascade, and 
raising the question as to whether the method could be 
extended to the entire blade for the cases of both pitch 
and plunge. Finally, Adamczyk [19], in a yet-unpublished 
work, employs the Wiener-Hopf method to obtain the pressure 
distribution for the linear cascade with subsonic leading 


edge. 


1.3 The Nonlinear Subsonic-Leading-Edge Cascade. Attempts 
to deal analytically or numerically with an oscillating 
cascade of blades having thickness are not nearly so numerous 
as their linear counterparts. Carrier [20] provided an exact 
theory for a single oscillating wedge. VanDyke 21,228 

also treated the isolated wedge as well as an oscillating 
КОШЕ ОТ Of arbitrary profile.  Teipel [23,24] SS cd an 
ingenious method-of-characteristics approach to numerically 
obtain the unsteady pressures on a single oscillating profile 
of arbitrary shape, obtaining good agreement with the single- 
blade solutions of Carrier and VanDyke.  Teipel's method 
assumes potential flow with negligible entropy changes 

across leading-edge and trailing-edge shocks.  Kurosaka and 
Edelfelt [4] have obtained numerical results for an isolated, 
oscillating, parabolic-arc airfoil and are presently working 
BEN ud the3r solution to a cascade of such airfoils. 
Chadwick [5] employed Teipel's method to obtain the unsteady 
peescmaeweistribution over the preinterference zone ОЁ ап 


arbitrary blade in a cascade of wedges. His results indicate 





that nonlinear thickness effects significantly alter the 


unsteady pressure profile of an arbitrary blade. 


J.4 The Present Investigation. The present work has two 
major parts: The first is analytical and treats the linear 
cascade while the second is numerical, dealing with a non- 
linear configuration. 

In the linear effort a closed-form analytical solution 
for the pressure-difference distribution along a blade in 
a Slowly-oscillating cascade is obtained. The solution 
technique is that developed by Chadwick [5] in the pre- 
interference zone (figure 1.2.1) of a pitching cascade, in 
which Sauer's solution of the unsteady linearized potential 
equation is employed. In the present work Chadwick's solu- 
tion is extended through all reflection zones and into the 
blade's wake. In addition, the solution is obtained for 
two degrees of freedom (pitch and plunge). Two standard 
configurations — Verdon's cascades A and B [12] — are con- 
sidered, one from a finite-cascade approach and the other 
as an infinite cascade in which Kurosaka's passage approach 
is employed. Extensive comparisons are made with existing 
solutions [6,12,19,17]; and several special topics, including 
cascade resonance and a unique inflow condition, are briefly 
treated. 

The nonlinear effort is an evaluation of the feasibility 
of extending Teipel's unsteady method of characteristics 


Eomemedseade of airfoils with arbitrary profiles.  Teipel 





K represents the nonlinear flowfield by a harmonic, 
unsteady, velocity potential superimposed on a steady one, 
and discontinuities in the flowfield variables across the 
shockwaves emanating from the airfoil are determined from 
the isentropic Rankine-Hugoniot relations. These concepts, 
coupled with the characteristic grid and compatibility relations 
derived by Teipel, permit determination of the steady and 
unsteady preinterference pressures acting on the second 
blade of a cascade of parabolic-arc airfoils. The numerical 
results are compared with experimental data from tests con- 
ducted with a pitching cascade by Fleeter and Riffel [25] 

and with linear theory. An improved method for determining 
the jump conditions across the bow shock is presented and 
the computational procedure is carefully outlined to facili- 
tate extension of the technique to the entire flowfield of 


a completely general cascade. 


1.5 Some Comments on Notation. In the linear portion of 
this investigation it is the unsteady, dimensionless, 
pressure-difference distribution along a particular cascade 
DieGemthacems sought. This pressure difference is a complex 
quantity given the symbol AC, and represents an amplitude 
function of distance along the blade. The harmonic oscilla- 
tion factor m is always implicitly understood. In 


addition it is desirable to obtain an expression for 


Le, at the instant that the blade in question achieves 





maximum up displacement, referred to as the initial position. 
This expression is obtained by multiplying through by the 
-і р ee — -ip 
factor e and so it is indicated by the sumbol Асе : 
In the nonlinear portion of this study only the local 
Bressures — Co upper and lower — are sought. They are 
also complex quantities and it is convenient to express 
them as a magnitude and phase angle relative to the motion 


of the blade under study; the symbols Kä апа Sp ; 


respectively, are employed. 
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2. THE LINEAR LOW-FREQUENCY CASCADE 


2.1 General. A cascade of flat plates with chord c 
oscillating harmonically with frequency & іп а steady, 
uniform flow U is said to oscillate at a dimensionless 


reduced frequency of oscillation k defined by 


n d EN 





er x? is small when compared with k , then the cascade is 


termed a low-frequency or slowly-oscillating cascade. Although 
low-frequency flutter is not common, it is worthy of investi- 
gation for several reasons, the major one being that it'is 

the only case for which one can hope to attain a relatively 
simple, closed-form expression for the unsteady pressures 

on a given blade. Such expressions provide insight into the 
flutter phenomenon, as evidenced by the identification of 

a unique inflow condition in a subsequent chapter of this 
Study. Also, Te low-frequency cascade provides an important 
limiting case for any general-frequency solution, numerical 


or analytical. 


2.2 Problem Formulation. 
2.2.1 The Linearized Potential Equation. The flow 
variables in the cascade flowfield are assumed to be 


composed of a freestream component and a much smaller 





perturbation component 


r-t +P pou 


The cascade is assumed immersed in a perfect gas, each 
particle of which maintains constant entropy as it moves 


through the flowfield 


DS 


DE =0 (2.242) 


In addition the flow is assumed inviscid and irrotational, 
guaranteeing the existence of a velocity potential 


WERT O2 


The continuity and momentum equations for the flow are, 


respectively 


“+p 7:7 =0 (ото) 


Dt 
IDEE 
B qs (2.2.5) 


Introducing the velocity potential and intesrating (2.2.5) 


produces, after linearization 


og 9% Р-Р. 
— — + = 3 232.36) 
St +U == Po ( 


Using the equation for the freestream speed of sound 


2 Ro 
= Po 29287. 
As > PE ( ) 


the continuity eauation can be linearized to give 





35 | 395 | —_ [2 + v 2—] =o (2.2.8) 
ad Ai 2 22 f 9s 9t. 9X 
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From (2.2.0) one obtains 


ә? _ o ә 
Se 7 f = +u Se) (2.229) 











and 


ә? ag ENZ 
ax pe ees FU Se] 22-100 


substitution of (2.2.9) and (2.2.10) into (2.2.8) then 





produces the linearized unsteady potential equation in 


two dimensions 











A M e$ ` a DD 
05941 a axat E IMD + ay? on 
where M is the freestream Mach number 
15 
МЛ == E ы (2.25125 


2.2.2 The Linearized Pressure Coefficient. A generalized 


form of Bernoulli's equation follows directly from the 


momentum equation (2.2.5) 
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= 2 P 
op KA ӨР _ 
5t d ч — CON ST. (2222129) 
foo 
Introducing the perturbation quantities in the form of 


C and ignoring the products of perturbation 


quantities one obtains 


ap ap Е ааа | 
SE + Us NS zo (2.2.14) 


where the prime notation has again been dropped. Thus, the 


dimensionless pressure coefficient is given by 


"Ra Fa 


Dur" 
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The Linearized Boundary Conditions. The are 
of a body placed in the flow can be described as a function 


of space and time 
F(N.+)= O (2.2.16) 
where r is the position vector of points on the surface. 


The condition that the flow at the body surface be everywhere 


tangent to the surface can be expressed by 


O 

T] 
i 
О 
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The plunging flat plate. The pitching flat plate. 





Figure 2.2.1 Two-degree-of-freedom oscillation. 


Berne surface. If the body is a flat plate oscillating 
harmonically in plunge and/or pitch, as shown in figure 2.2.1, 


then the describing functions are respectively 
Et eR 5 (2.2.18) 
Eg Oz yr rn eto сог, 
where h. and Eë are the maximum displacements and b 
is the pitch axis location on the x axis. Applying the 


generalized flow tangency condition (222.] 7] СО СОЖА 18) 


and (2.2.19) yields the boundary conditions for plunge and 





pitch respecitvely 


ə 
2 


=—hiwet** at PLI (2220) 


U 


= o U + O ^t otyz=o (2.2.21) 
AA Nond»menstonalization. It is convenient to intro- 

duce nondimensional variables into the governing equations 

to simplify and generalize them. The distance variables 

(x,y,b) are normalized by chordlength c , velocity potential 

Br che product Uc and time by the quantity c/U . It is 


also convenient to set 
„>= © Xas L RADAN ` (052.22) 


and to separate the velocity potential into a complex 


amplitude and phase angle 
LRE 
(x o t = pape 22,908) 


Thus the potential equation (2.2.11), pressure coefficient 


(2.2.15) and boundary conditions (2.2.20,21) become 





19 L Ë рм + 2 w 4222 oO 2282 
P >x da" к = | 


— -2 (2b. DDR 
= RD) ( ) 





s 


39 A 
2D 
ag = Tl1+ iria] at y=o (2.2.27) 


where, 


ZS: M*- 4 (2.2.28) 


апа es and Ф аге complex amplitude functions of x and 
y only with the factor git omitted. 

ARS Cascade Geometry. The physical parameters chosen 
to describe the cascade are the Жы л теб T between 
the blades d and the horizontal distance between the 
leading-edge Mach waves m. Both "m are normalized 
by chordlength. The cascade can also be described by 
solidity Z and stagger angle O6 , which are shown in 


figure 2.2.2. These parameters are related as follows 


mue: sis (2. 2.29) 
SIN 9 

= — —— 2.2.30 

$ 7 ( 30) 


where А is given by (2.2.28). The linear cascades employed 
in this study are Verdon's cascades A and B [12], which 
have become standards in the field and are depicted in 


ште 2.2.2. 





Zao The Isolated Airfoil. Consider a single flat plate 
oscillating with two degrees of freedom as shown in figure 


2.2.1, and define 


S-X-87 Djs (25501) 


Then s is constant along any leftrunning (upper surface) 
Mach wave, while s is constant along right-running (lower 
surface) waves. Sauer [18] showed that the potential 
equation (2.2.24) has the following left-running and 
right-running solutions, respectively, for a slowly- 


oscillating airfoil 
о . UM? 
Pg guia + e [nus - une (203802) 
" ту, с” 2 

Dr (%, 9) = 9 Sr JN ER A 9.03) (293.3) 
where the g and h functions are determined from the 
Beundary conditions (2.2.26,2/). Requiring that the 
velocity potential be continuous across the leading-edge 
Mach waves, one has 


glo) = noO = ©.(0> = (9) = О (2.3.4) 


Enforcing flow tangency at the upper surface of the airfoil 


yields, for plunge and pitch respectively 
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PL: Jat =O б о е (2.3.5) 
А 
PI: (== biens Si bas) (2.3.6 
g^ A = E А? 2 P (2.3. ) 
And so on the upper surface 
PL: drun, ais SÉ s (2 37) 
2 
x RTL CI g? 2 a 
ET: Gi (%,9) А ALS + step + вм) 238) 


Thus the upper-surface pressure coefficient follows from 


(2.2.25) 
PL: š сё = – 66 (2.3.9) 
š o _ _ 2 2R ES 2 
PI: 5 5 |5 supra] (2.3.10) 


where the subscript indicates the single airfoil. Similarly, 
for the lower surface 

> 

St 


=~ Cp (2.3.11) 


2.4 Conermurcy of the Neto sie Potential Across a Mach Wave. 


Consider two reflection zones separated by the Mach wave 
s = So , as shown in figure 2.4.1, and let the variables z 
and Z respectively measure distance tangent and normal to 


the Mach wave. The velocity potential in zone II is simply 


As 












“Mach wave: SzSg 





FLat-piate airroil 


igure 2.4.1 Variation of the velocity potential along 
the boundaries of reflection zones. 


at in zone I plus an additional disturbance 
Q" = ét + YP (2.4.1) 
nce the tangential velocity component remains unchanged 


ross the wave, (2.4.1) implies that Pa 15 а Sunctilon oF 


alone, along the Mach wave 
a= HZ) For S= So (2.4.2) 


t since Z is itself constant along the Mach wave, then 


| must have the form 


p= CONST. for S=Se (2.405) 





That the constant is zero follows from the physical con- 
straint that the velocities remain finite in the wave 


meself, and so 


фс. = фЕ (ъл ОАЕ) 


which is in agreement with the conclusion reached by 


Landahl [26]. 


2,5 The Finite Two-Degree-of-Freedom Cascade. As an 


elementary means of introducing the method of analysis, 

the pressure-difference distribution is found for the second 
blade of a finite version of cascade A. This analysis also 
affords the opportunity to compare the finite-cascade results 
obtained with numerical finite-cascade data [17] and SE 
tical results [6] of existing solutions. It is emphasized 
mat the pitch results through equation (2.5.27) were first 
obtained by Chadwick [5]. The reflection zones referred to 
weeen e tured figure 2.5.1, and the pressure results are 
given in terms of the isolated-blade pressure coefficients 
(2.3.9,10). Before proceeding with the analysis, some 


additional variables are introduced for blade two 


v= $- A? S= +p (2.51) 


which are related to those of blađe one by the translation 


relations 





Cascade A (not to scale) 


_blade 3 
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Figure 2.5.1 The finite-cascade reflection zones. 
X = € + m. д у= 9+9 (2.5.2) 
ъ= ҝет 5 = 6+ п + 244 (2.5.3) 


2.5.1 Zones I and IA. In zones I and IA Sauer's single- 


blade solution applies 


DT Lui, 
Фес уу = 9 +e[lno- gg] (2.5.4) 


4 4 
DEJA Rira- T| (2.5.5) 


and in zone TA the boundary conditions take the form 





IA 
PL: 39: 








Sq Ek e at %=0 (2256) 
ә@ї^ | 
ET: SCH = -|1 + ik -h| e` at 2=0 ЛБ 


where ps is the interblade phase angle by which the second 
blade leads the first. Applying the boundary condition at 


the upper surface of each blade and using that 


ON = nod = 0 (25sec) 


gives, for plunge and pitch 


PL: 905 = О Ф) = А (2.5.9) 

Е = улс» 2 5510 
PI: fie) = E 0 52 е“ + р (2.52.10) 
PL: nS = 5 nas is e +С (2.5.11) 


: 2 
hs es UE Cr ets) 
PI: (255: 12) 
h. (97 = — GE pato) ein ерам" +D 


miere A, B, C and D are constants of integration. By 
enforcing continuity of the velocity potential across the 
leading-edge Mach wave of blade two, the constants are 

evaluated and the velocity potentials determined in terms 


of the isolated blade potential Фү (2.3.7,8) 





PL: pr = fr I^. prev + rds (es ey 


Фе = Ф 
PI: (255 4) 
EL A Gua ac 


where £ is a constant that need not be evaluated because 
it becomes the coefficient of a term in Re when the 
pressure coefficients are Er, Thus, on the upper 
surface of the second blade, in terms of the isolated-blade 


pressure coefficient Се, (22350718) 


PL: en се” Pal) 
PI: AN er + Em (2.5.16) 


2.5.2 Zones IIB and VIB. In zone IIB the velocity 


"potential is that of the isolated blade (2.3.3) 


p= D EN 


while in VIB the potential acquires an additional 





Ihe symbol @ is used throughout section 2 to indicate 
a constant which need not be determined. 
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rieht-runnäng: disturbance 
"s 
Bae = qz + 96(S) + | he (5) + gta! (2.5.18) 
ps3 Zone Il. In zone II a right-running reflection 


occurs and, using the translation relations (2.5.2,3) to 


shift to the blade-two coordinate system, one has 


ES : = ¿mi = 

PL: ФЕ = Biv gh) + [= + h, (3) + glo 2.5219) 
R | 

бз = + = = М, ЕГА 

BT: (225.20) 
SSC > «|n, Gà e = 2 2.0] 
where 

sto) =h,(o) =O (2253.21) 


Applying the flow tangency condition at the lower surface 
of blade two and using the Sauer functions to construct the 


velocity potentials yields 


PL: De = gi een em] (295022) 
= = i SENE em 
E EE ж С ^у— jm (V - &*^) 
PI: Slam + Рам? + Co g*- gv? py - ev) (2.5.23) 


m аА л CG? Es 





And the pressure coefficients follow from the potentials 


PL: Се, = Cp 12-e”) rs) 
РЇ: СЕ cg, 2- €) 555 Un- zj'm-284M) (2.5.25) 


Thus the pressure difference can be computed through zone 
E and by multiplying by the factor e "/* the pressure- 
difference amplitude is obtained at the moment the second 


blade reaches maximum displacement 
PL: Al, = -2 co G — е“) of € <im (225626) 


AC e- = -208 0- m pam A "m) e 


PI: (2.5.27) 
© <\-т 


one TITI. In zone 111 a left-running wake distur- 
bance is experienced, and after applying the translation 


relations 


m? 
PL: De = 9? er E yg | 2292208) 


Dez DE — ae = S iE cm- 28d-- gdM*)0 - ev) 
p EN es „(орам рима es] (225,29) 


-P et is (S) + Ww SS ads e| 





where 


9:0) = зб) = О (55550) 


Since zones III and VIB are separated by the wake of the 
blade rather than the blade itself, the wake boundary con- 
ditions are used between the two zones. The first of these 


demands pressure equilibrium across the wake 


u CF at 9-0 (o s sn) 

which leads to the equations 
JOE = ST (Cp) No») 
ha cn - date t [ 549 90m] = Fm (CP) were) 


where the right-hand sides of (2.5.32,33) refer to the 


steady and frequency terms respectively of e (2. 329510) 


р1 


The second wake condition requires continuity of the normal 


velocity component through the wake 


nr VIB 
әй _ 927 at 9-0 (2.5.34) 
92 od 


which leads to the equations 


^(S) + ek (S) = O (255 35) 
dst 





: z 
hy (5) +h CS) + Gi E =© EE EE ER 


mations (2.5.32,33,35,36) are solved for the zone-III 


disturbances 
G: = ; acs 
P A) O h, (0) 2 (205557) 
° E eme em E = = 
p KK 2 (S—1) = Cy a (225458) 


2.5.5 Zone V. In zone V an unknown disturbance is added 


o? = Pr +93(0) + 9s CS) + к} һу + hs (©) 
(2.5009) 


LM" 


z ( EX vo — ge 


Bo 


which is determined by enforcing flow tangency at the lower 


surface of blade two in zone V 


d 


PL: As) =° ns() =~ 4 (2.5.40) 
I = E 
РТ: (2.5.41) 
hs (©) = e BT "kg éd M?) 


which leads directly to the pressure difference for the 


last portion of the blade 





PL: AC, e^ z-20p i-m«€ € I (2.5.42) 


* Y * ` 
ACp eK = cn ECH, e A 
PI: (2.5.43) 
l-m <§ 21 
2.5.6 Results and Conclusions. The finite results for 


the second blade are compared with four existing solutions 

in figures 2.5.2 and 2.5.3 for three interblade phase angles: 
Kurosaka's solution [6] is a closed-form, analytical, 

ce cescade solution; Platzer's [16,17] is a numerical, 
 method-of-characteristics solution for the finite cascade; 
and Adamczyk's [19] and Verdon's [12] solutions are numeri- 
cal, infinite-cascade results. It is interesting to note 
that the Free en: agree exactly with Kurosaka's | 
infinite solution over the entire blade for the plunge case 
and through zone II for the pitch case. Over the zone-V 
portion of the pitching blade — the portion exposed to the 
wake disturbances from blade one — the two solutions agree 
numerically, as seen in figure 2.5.3, but are not analytically 
the same. In fact, Kurosaka's solution is singular for 

Kr О (in>=phase oscillation of the blades) over the last 
portion of the blade, while equation (2.5.43) behaves well 
for that case. All five solutions were found to be in 
excellent agreement for interblade phase angles larger than 
15°, but for small angles considerable disagreement is 


evident, as seen in figure 2.5.2. 
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Figure 2.5.2 Comparison of the finite-cascade 
solution with existing solutions. 
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Figure 2.5.3 Comparison of the finite-cascade 
solutions with existing solutions. 





Chadwick [5] compared his finite and infinite results 
for the pitch case over the in-flow portion of the blade and 
noted that although the finite and infinite pressure- 
difference Ас) distributions are the same, the local 
pressure distributions (Co: с) may in fact differ. The 
author found this phenomenon occuring over the remaining 
portion of the pitching blade also. Figure 2.5.4 shows that 
the local pressure distributions for the present finite 
solution and Verdon's infinite solution do not agree as 
well as the pressure-difference distributions. Some con- 
clusions can be drawn from the analysis at this point: 

1. For low frequencies of oscillation the finite cascade 
rapidly approaches the infinite-cascade solution, achieving 
good Me ent with the pressure-difference distribution 
of the latter at the second blade. 

2. For small interblade phase angles, the present solution 
agrees with Kurosaka's Laplace transform solution for the 
infinite cascade but considerable differences with the other 
published solutions are noted. 

3. Although the present finite-cascade solution and 
Verdon's infinite-cascade solution for the pressure-difference 
distributions are in agreement for large interblade phase 


angles, the local pressure distributions resulting from the 


two solutions are seen to disagree. 


2.6 The Infinite No Degree „or sEreedom Cascade. In order 


to demonstrate the applicability of the above analysis to 
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Figure 2.5.4- Comparison of the finite and infinite 
local pressure distributions. 


the infinite cascade, it is applied to an infinite version 
Of Cascade B. Since the advantage of an identifiable first 
blade is lost, a condition specifying the space periodicity 
of the flowfields between the blades is derived. The pitch 
results through equation (2.6.17) were first obtained by 
Chadwick [5], and the reflection zones referred to in the 


solution are depicted in figure 2.6.1. 
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Figure 2.6.1 The infinite-cascade reflection zones. 


2.6.1 The Flow Periodicity Condition. The periodicity 


£ the flow from blade passage to blade passage was first 
dentified by Kurosaka [6] and is the heart of his passage 
pproach to the infinite cascade. An alternate derivation 
ollows. 

Consider the corresponding zones I and IA, and assume 
hat the potentials are the same there except for a complex 


actor K 


$e OS gy = (eQ) ОШО 


ince each potential must satisfy the boundary condition of 


ts corresponding blade 














= =-iR at y=o c 
E = — k e SE n =o (2.6.3) 
fe is clear that 
K= e-m (26.4) 
апа so 
бл ох) = фал er ne 


The same result follows for any pair of corresponding zones 


on adjacent blades — for pitch 


ees low periodicity condition. 


and plunge — and is termed 


2.6.2. Zones I and IA. The analysis in zones I and IA 
proceeds much as in the finite case (section 2.5.1). The 
velocity potentials have the form 

T ¿Mm? ( 
ФЕ = 9.09) + к "GT Td 5) (2.6.6) 
2 
Pro = Jats) +R 2 u LA | (2.6.7) 
ß 
and the boundary conditions yield 
(2.6.8) 


PL: I (5) = А САС) = = 





PI: nis re JD = e+) (2.6.9) 


8 
° M? 
h(S) = + ° - y AS +E 
PL: (25610) 
nia ln = =o e tM? gerar 
A? 
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ns) == (E + bets + AMCs) + 6 
2? 2 
РІ: (2-6, ЕГ) 
nA(G)z "= [= == epo) e ^4 + mino] +H 


The constants A through D are evaluated using, respectively, 
continuity of the velocity potential across the leading-edge 
Mach wave of the (n+l)st blade and the flow periodicity 


condition 


@ (S=w) = @,, (=O) (256112) 
м=р е^ (2.6.13) 


The pressure coefficients, in terms of those of the isolated 


blade (2.3.9,10), follow 
PL: a (226214) 


g отт M v | 2«R m 5 
PI: A EP CL e) (2.6.15) 





m comes TI through IV. Finding the additional 
reflections through zone IV is accomplished by routine use 
Eu slow tangency condition (2.2.26,27), and so only the 


pressure results for the (n+1)st blade are given below 


PL: сре‘ = -2 сро- ет“) O£4 <2Ad (2.6.16) 
ACp €^ 2-20, Q- ev) FR (gam ptm” 
РЇ: deis jen) 
GER: — 


PL: Ateg = —4- me" s, kom (2.6.18) 


SEH (2.6.19) 


E k )- m 


ЭИ ООО е5 у, Vr, and VII. In the first three zones 


above the wake of the nth blade 


Di = da pas ВЭБ R Insta - yg! (26.20) 
Xy 
DZ- AFX agr nat S 1940) ОО о; 
К+ $ 
EN SS Seel [nion fra) (2.6.22) 


where the 94 апа h Sauer functions are known from the 


= 


analysis in zone IV and 





Jan ns() =O (2:00:23) 


9; C- mz hz 1-и) = О (2.6.24) 


2.6.5 Zones VA and VIII. In the zones straddling the 


wake of the (n+l)st blade, using flow periodicity and 


results from zones IV and V 


ZA DIA 


eo 
nei = Poa + Isa (NM + heat) t дви (2562525) 


+ кўп? es + hy (07) + hen) 


LM? 


Lo sd 


(2:56:25) 
N : 2 ` ` | 
A Е е^ 9 G)- 2] _ 54 gm = 


where 


aa 
gar) = «| nsa teo = `> i 3st) = 


E | 
190 + R ЕЕ оре so) t ew 


(2.6.27) 


930 = n (=o (2.6.28) 


The wake boundary conditions — pressure equilibrium and 


normal velocity continuity — are enforced petween zones 


VA and VIII 











NA _ zu 
Се = Cr, at t=O (2.6.29) 


YA 


ə Spo 
~= = at {шо (2.6.30) 


giving rise to two equations in three unknown Sauer 


functions 





9su-e-)-9 — 9, = ST E ЕР 5 сонлу (2.6.31) 
920-67) + arg =0 (2.6.32) 
hs li~ €777) - V, ие + —- d 95 QE 
| (2.6.33) 
ei [sse en e =н A 
hs(1- Se, + SS E Late e 2o ger págs enn] =° (2.6.34) 


where the right-hand sides of (2.6.31,33) indicate the 
steady and frequency terms, respectively, of the enclosed 
quantity. The third equation in the set comes from applying 


НЕ 


the boundary condition to Gë at the lower surface of 


the (n+1)st blade 
Js S -gy =O (2.6.35) 


hs e7- M at Ys Qe 7 — 89 95 еги) =0 9.6. 36) 





Solving equations (2.6.31) through (2.6.36) for the Sauer 
functions permits determination of the pressure difference 


Meer the last portion of the blade 


PL: Sener 2x6 im < « 28d et сЕ) 
As er nen ea n 
ES (2.6.38) 
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PL: Ate ECH, Lë El 84+ т< { & \ (2.6.39) 


Асрет = -2 Ci (2- e^) — 512 (pdm 22 m) 





PI: + [apm Е) Q- e`) re, (2.6.40) 
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2.6.6 Results and Conclusions. Figure 2.6.2 shows that, 
as with the finite-cascade solution, there is considerable 
disagreement among the several existing solutions for small 
interblade phase angles, but even for those angles the 
present solution and Kurosaka's are in agreement. The 
pressure-difference expressions derived above for the infinite 
cascade agree exactly with those of Kurosaka [6] over the 
entire plunging blade and up to £= 2pd+ m for the pitching 
one. Over the last portion of the pitching blade the two 


Solutions are slightly different analytically but are in 
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Figure 2.6.2 Comparison of the infinite-cascade 
solution with existing solutions, 
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numerical agreement as seen in figures 2.6.3 and 2.6.4; 
the analytical discrepancy could not be resolved. All of 
the solutions are in accord for the larger phase angles. 
Continuing with the conclusions begun in section 2.5.6, it 
Can be added that: 

4. For interblade phase angles larger than 15? the 
present method of analysis yields closed-form pressure- 
difference expressions which are in agreement with existing 


solutions for the slowly-oscillating cascade. 


2.7 Limiting Cases of the Finite and Infinite-Cascade 
Solutions. 


2.7.1 The Steady Cascade. А11 terms of the plunge solu- 
tions are frequency terms; therefore if k is set to zero, 
the Ee enee distribution goes to zero for She 
entire blade, which is, of course, what should happen for a 
stationary flat plate aligned with the flow. 

If k jis set to zero in the pitching case, then one 
obtains for the isolated airfoil 


Cp, =- (2. 2:905) 


“м 


and so both the finite (2.5.27) and infinite (2.6.17) 


solutions becomes 
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Figure 2.6.3 Comparison of the infinite-cascade 
solution with existing solutions. 
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Figure 2.6.4 Comparison of the infinite-cascade 
solution with existing solutions. 
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where the phase-angle term has been ignored because it is 
meaningless in the steady case. Equation (Raid ts. the 
classical Ackeret solution [27] for a flat plate at an angle 
of attack equal to one radian, which is the maximum up-pitch 


Position. 


2.7.2 The Sinqle-Blade "Cascade". If either cascade is 


"squeezed" down to a single blade by setting the spacing 
parameters m and d to zero, then the solutions (2.5.26,27) 


and (2.6.16,17) reduce to those of the isolated blade 
ACp= -2.Cp (24470) 


where the phase term has again been ignored.. 


2.7.3 The Infinitely-Spaced Cascade. The isolated-blade 


solution is again recovered from equations (2.5.27) and 
(2.6.17) as the cascade spacing is allowed to get large. 
If the spacing parameters m and d approach infinity in 
such a way that the terms pam? and рап remain equal, 
then equation (2.7.3) again results. The significance of 


these two terms is discussed further in the next section. 


ENEE Inflow Condition. A consideration of the 


1 
infinitely-spaced cascade in section 2.7.3 leads one to 


1тһе su orris grateful to Professor M. F. Platzer for 
mn s n Ing this question. 





«cz For which cascades does the relationship 


(2. 75» 


aM’ = gum 


om equation (2.5.27), for such a cascade 


11d? As seen fr 
on of the blade 


te pressure-difference over the inflow porti 


i merely that of the isolated pitching blade times Ene 


Bctor 0- e. 


Let & be the Mach angle defined by 


X= SINT a (Qe) 


Men it is not difficult to show that for any cascade having 


he stagger angle 


acer су л (2.7.6) 





iquation (2.7.4) “n also hold: Saia another way, all 


nembers of an infinite family of pitching cascades with the 
stagger angle given by (2.7.6) — and with no restrictions 


5n blade spacing — will have exactly the same inflow pressure 


distribution, i.e. 


| 


PI: Acep E %) Een 


This is termed, by the author, the unique inflow condition 


¡for the slowly-oscillating cascade. Figure 2.8.1 depicts 
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Figure 2.8.1 The unique inflow condition. 





three different cascades, each with the same stagger angle, 
satisfying (2.7.6). The pressure distribution for each 
was obtained via Platzer's method-of-characteristics pro- 
gram. The unique inflow condition is evidenced by the 
exact agreement of the pressures for all three cascades 


over the first portion of the blade. 


2.9 Some Comments on Resonance. For a given reference 

blade in an infinite cascade, resonance occurs when aero- 
dynamic disturbance waves generated by upstream blades are 

in phase with the motion of the reference blade [12]. 
Samoylovich [28] shows that at the resonance condition the 
linear cascade experiences zero anio damping. The 
Mevelopment that follows is based on Verdon's [12] excellent 
explanation of the resonance phenomenon and provides a means 
of delineating what he calls the subresonant and superresonant 
regions of the a plane. 

Referring to figure 2.9.1, also taken from [12], it is 
seen that the leading edge of each blade creates disturbances 
which propagate in all directions at the freestream speed 
of sound аы. These disturbances are swept downstream by 
the freestream flow U , and there are exactly two resultant 
waves which propagate along the cascade leading edge in 
the direction opposite to the direction of rotation and at 


velocities Vi and У. . From the law of cosines one has 


i 
М M(m+84)17 [m(m+ 21] “ (2.9.1) 
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Figure 2.9.1 The backward-traveling waves 
causing resonance. 
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Figure 2.9.2 The subresonant and superresonant regions. 





The times tin! ton required for a wave from an upstream 


NEE E, -2, ...) to reach any given reference blade 


(п = 0) are 
-nM % 
tian = в? Mon a) [min 24a] ? (2. 9009 
The reference blade lags an upstream blade by 
n 
Ka "e nz--2,... (28503) 


and resonance occurs when the upstream waves are in phase 


with the reference-blade motion, i.e. when 
o7 bin . (2.9:4) 


or when 





AT LIE то} (22943) 
For a fixed reduced frequency k the superresonant inter- 
blade phase angles lie between J^ and Mos as shown in 
figure 2.9.2, while the subresonant phase angles lie to 
either side. The superresonant region is important because 
Verdon is of the opinion that, in that region, neither a 
finite-cascade solution nor a passage-approach infinite- 
cascade solution will adequately approximate the true 


infinite-cascade solution. In an effort to shed some light 





on the subject, some comparisons are presented in figures 
EN cthrough 2.9.6. For cascade A, with k = 0.05 , the 
superresonant range of phase angles is approximately -7° 

to -2° „ and so -5° was examined. Figures 2.9.3 and 2.9.4 
show that Platzer's finite solution and the author's passage- 
approach infinite solution are in agreement with Verdon's 
solution. For cascade B, however, Verdon's solution is 
significantly different from the other two over the last 
portion of the blade for the superresonant phase angle -3° . 
Figures 2.9.5 and 2.9.6 show some high-frequency comparisons 
for the phase angle -90° , which is in the superresonant 
regions of both cascades A and B. In these two cases the 


two solutions — one finite and one infinite — are in accord. 
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pressure distributions. 
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Figure 2.9.4 Low-frequency superresonant 
pressure distributions. 
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Figure 2.9.5 High-frequency superresonant 
pressure distributions. 
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3. THE NONLINEAR CASCADE 


3.1 General. In two papers Teipel provides the numerical 
tools for the analysis of low supersonic flow past a single 
Es sting airfoil of finite thickness. In one paper [23] 
he presents an unsteady, nonlinear characteristics method 

to determine the continuous variation of flow variables. 

In the other [24] he treats the oscillating shock waves 
present in the nonlinear problem to permit determination 

of the discontinuous jumps in the flow variables across the 
shocks emanating from the single blade. Using these tech- 
niques as well as some computational and theoretical insight 
provided by. Chadwick [5], who very ably TE and applied 
Teipel's work, the present study is aimed at finding the 
pressures, steady and unsteady, on the second blade of a 
parabolic-arc cascade. The cascade oscillates in pitch 
alone, and the pressures are sought over the entire upper 
blade surface and over the lower surface up to the point of 
impingement of the trailing-edge shock from the preceding 
blade (figure 3.2.2). These surface pressures are computed 
by means of a FORTRAN IV computer program prepared by the 
author for the IBM-360 computer. The program listing appears 
at the end of this work immediately after Appendix A. What 
follows is an orderly presentation of the theory underlying 
the computer program and the computational techniques used 
in implementing that theory. Wherever possible the text 


notation will coincide with that used in the program. 





Baez) Cascade Geometry. 

3.2.1 Blade Shape. The blade shape was chosen to approx- 
imate the blades in the test cascade employed by Fleeter 
and Riffel [25] who attempted to measure steady and unsteady 
pressures in a pitching cascade. The airfoil used in the 
computer model has a flat upper surface and a lower surface 


conforming to the parabolic arc given by 
ЕЕ XUX=1) (332230) 


where chordwise distance x and midchord thickness "c are 
normalized by chordlength. Figure 3.2.1 shows the computer 
model for t = 0.03 superimposed on Fleeter's test blade 
and lists the surface coardinates of the test. blade in per- : 
cent chord. The greatest discrepancy between the two shapes 
occurs on the aft portion of the blade in the region behind 
the lower-surface shock reflection; however, the lower- 


surface pressures are not computed in this region. 


Е 2.2 ризасшэрасапа The blade spacing parameters are 
input variables for the program and are discussed in more 
detail in section 3.5.1. The physical spacing of the blades 
of the computer-model cascade — dubbed Cascade T — is shown 


moe gure 3.2.2. 


3.3 Teipel's Method of Characteristics. The underlying 


insight of Teipel's treatment of the unsteady problem is 









Test bDlade shape — — —— — 
Computer model shape 
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` 0.00 т тоз чш — — 
| UPPER SURFACE ` | LOWER SURFACE 
г р 
0.00 -0.26| 49.355 0.00|| 0.00 -0.26| 50.51 -2.97 
5.21 -0.14| 55.65 0.10|| 6.29 -1.08| 56.83 -2.84 
Wr s 61.95 °0.191||112.59 -1.65| 63.15 -2.60 
EEN 16 68.25 0.291118.90 -2.12| 69.47 -2.27 
26.14 -0.15] 74.55 O0.411]25.21. -2.49| 75.78 -1,83 
30.45 -0.13| 80.88 0.50||31.53 -2.76| 82.06 -1.37 
36.755 -0.10| 87.21 0.26||57.85 -2.95| 88.36 -0.96 
43.05  -0.04| 93.53 0.26||64.18 -3.00| 94.66 -0.62 
99.85 -0.10 99.85 -0.10 
* Expressed in percent chord. 
Surface coordinates. j 


Figure 3.2.1 Test blade shape and coordinatess 
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Figure 3.2.2 Geometry of the computer model. 


i representation of the velocity potential as the sum 


fa steady component and a small-amplitude unsteady one ` 


PILYE = PAg dou) e * (ЕСЕ) 


depicted in figure 3.3.1. When (3.3.1) is substituted 


ko the transonic unsteady small-disturbance equation 





p S6 ro mob Mito (3.3.2) 


a 1 or ME 
[ F Er SX SX: EE эх ә a+? 


id into the unsteady boundary condition 


lror a discussion of the limitations imposed by this 
luation see Teipel [23] апа Chadwick [5] 
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Figure 3.3.2 The characteristic grid. 
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where y(x,t) describes the oscillating airfoil surface, 
one obtains two separate problems coupled by the steady 


velocity potential 22 . One problem is nonlinear, steady 


and zero-angle-of-attack 





9% _ 92 


and the other linear and unsteady 


ai Ié Dy 39. 5% 219% az | 
[e MS ee 2 + mas ` + 2¿ uff ae б фе С] 


od. 


E ie cox e| OSET) 


where D (X,Y) is an amplitude function and the factor Qt 
is omitted. It is important to note that since the char- 
acteristic directions are determined by the coefficients of 
the highest derivatives, both problems will have the same 


Т 
characteristic directions. The state variables” for the two 





Ihe symbol 44 is used both for the steady state variable 
and the interblade phase angle, and it is normally clear 
from usage which meaning applies. 





problems are essentially the steady and unsteady velocity 


components u. , v , u and v 


О О 1 1 
STEADY: A= Pr +M*(S-41) Ue (3.3.8) 
B= M2 MES +1) Vo (3.329) 
UNSTEADY: Ur, Ya (3.3.10) 
where 
fM (seo) 


Encochbaracteristic directions for both problems are: 


QO 


s< 


Mod 
SL? o б, 





where the upper sign refers to the leftrunning & charac- 
teristics and the lower to the right-running A character- 
ШО сз." The compatibility relations, however, are not 


the same for the two problems. They are 


STEADY: Mae Ka CONST, в (ER) 


2 ,.4 
Lohe symbol 8 is used both for the quantity (M -1)* 
See Trigt- running characteristic direction. It is 
normally clear from usage which meaning applies. 











du, T N ES Bap + x i dx 
UNSTEADY : мг (373 11) 
M 
O 
where 
Ay=22) ax Б = dx (3.3.15) 


P 


Equation (3.3.14) requires some additional comment. 
Note first that in addition to the state variables ч, 


апа Vis g4 is also present in the unsteady compatibility 


relations. Since 


SÉ Di y . ; 
ad, = Agr A == y dy (3.3.16) 


then Йй, must be found by integration as one moves from 
one gridpoint to another. For the points A, B, and C 


shown in figure 3.3.2, for example 
с 
ax: фо = %“ + Js — 4x (3.3.17) 
А N * 
S Mi 
8 —— M 
А: = 9 + Jo р (333.298) 


In performing the actual computations an average of the two 


above expressions is employed. 





The factor 2 в SS AS straightforward, repre- 
senting the increment in ^ along the А characteristic: 
At point C in figure 3.3.2 for example one could use the 


approximation 


Ag = Ar ^ь (3.3.19) 


However Ax 1s a bit more complex, for while the increment 
dx is along the & characteristic, the quantity $x a 
represents the change in A with respect to x along the 


crossing Р characteristic at the point of interest. 


Again at point C one could employ 


>. № Ав te, | 
Ач = (X. — XA) (3.3.20) 


3.4 Teipel's Treatment of the Oscillating Shock Wave with 
Extension to the Present Cascade. 


3.4.1 The Isentropic Rankine-Hugoniot Pressure Relation. 
Consider a normal shock moving through an undisturbed 
perfect gas at velocity ve relative to the gas as shown 
in figure 3.4.1, where the hat notation indicates conditions 
behind the shock. The governing continuity, momentum and 


energy equations are 


except AO u. С and V., the velocities appearing 
in sections 3.4.1 inroügh DM are dimensional, in a brief 
departure from the normalization of all velocities by the 
freestream velocity U. 








Figure 3.4.1 Jump conditions across a moving, 
normal shock. 
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Figure 3.4.2 Consequences of the isentropic 
assumption. 





A A 
PENS = B Nn = Vs = Vin) (3 20900 


^ A 
P+ pVs = РВ Ф (У, м. Voy (2-402) 
VŽ MN y ES Өй 2 
n+ = E SEN) i Vn) (3.4.3) 


Assuming constant specific heats, the speeds of sound are 


introduced, and the enthalpies are eliminated by 








P 2. ^ 
o e H en д. = б-х- (3.4.4) 
Р 2 
ot ^ RŽ 
= n = 372455 
n In WS ( ) 


and, after a good deal or manipulation, the change In 
the normal velocities across the shock is found in terms 


of the shock velocity and sound speed ahead of the shock 





>` 
= О. 
М-М _ 2. = (RANG 





m (3.4.7) 


yield the Rankine-Hugoniot pressure relation sought 


P 2 у || Vs V 
ir LE) —1 (3.4.8) 





Figure 3.4.2 is based on figure 8-3 of [29], and it shows 
that the isentropic approximation introduces no significant 


error for Fleeter's test cascade or its computer model. 


3.4.2 The Blade-Two Shock Position. The leading-edge 
shock from the pitching blade oscillates out of phase with 
the blade. Let бу) be the shock angle of the steady, 
leading-edge shock as shown in figure 3.4.3. Since the 
shock is curved, e varies with distance y from the 
blade, and by convention 


ur 


OL Vo < (3.4.9) 


The shock angle at any point on the oscillating wave can 


be written as a small variation from the steady one” 
aint 
í = "Sek SN ES (3.4.10) 


In order to locate the oscillating shock relative to the 
steady one, Teipel introduces the complex displacement 
function g(y) , also shown in figure 3.4.3, and so the x 


coordinate of a shock point is approximately 


х= хо 2) COT Y. - a“ (E 





1тһе symbol gamma used here and in figures 3.4.4 and 
3.4.5 represents the total shock angle. All other uses of 
the unsubscripted, unsuperscripted yamma stand for the 
dimensionless ratio of specific heats. 
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Figure 3.4.3 The oscillating leading-edge shock wave. 
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Figure 3.4.4 Velocity components and the oscillating 
shock wave. 
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where (Xo Y) are the steady coordinates of the blade's 
leading edge. Taking the total differential of (3.4.11) 


and noting that along the shock 


әх 


ано ч (3.4.12) 
results in 
zl euni a, Së ett (3.4.13) 


dg 
thus relating shock angle and shock displacement. 


3.4.3 The Velocity Components Upstream of the Blade-Two 


Shock Wave. Ahead of the shock the velocity components 


have both steady and unsteady parts 
Um cuu e (3.4.14) 
vk = We went (3.4.15) 


where the barred velocities are dimensional counterparts of 
the dimensionless quantities uo, Vor U and Vi already 
mieeeduced in equations (3.3.8,9,10).. However, the flat upper 


surfaces of the blades of cascade T are aligned with the flow 


in the steady position, and so, upstream of the shock 


Uo = Yo =O (3.4.16) 





Referring to figure 3.4.4 one has for the normal апа 


tangential velocity components ahead of the shock 
Vn = U SIN Yo + US COS So S + (U SIN So + V, COS go) &'** (Saa 


V= U cos So - UI SINS 8? + (TH, COS Yo Y, SIN Yo) eint (3.4.18) 


3.4.4 The Blade-Two Shock-Motion-Induced Velocity. 
Relative to the gas ahead of it the shock is moving at 


the velocity 


Vg = —- Va + N! (3.4.19) 


where V' is a normal velocity induced by the motion of 


the shock, and is seen from figure 3.4.5 to be 


ı_ 9% 
MES SAN S (3.4.20) 


or, from (3.4.11) 


V'- QR SINS 9 et (3.4720 


3.4.5 The Velocity Jump Conditions Across the Blade-Two 


Shock. Using equation (3.4.13) to eliminate w' from the 
expressions for the upstream velocity components (340 1758). 


they become 
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Figure 3.4.5 The shock-motion-induced velocity. 


Vaz U SIN Wo + (LI CoS So wis ger SIN Yo + Vj COS So) ent (3.4.22) 





V.S U coss, «CU sit GS + U, coSs Yo = SINY) SE (3.4.23) 


Note that each of the above velocities has three component 
parts. The first term in each expression arises from the 
steady shock position, the second arises from the oscillatory 
motion of the shock about that steady position, and the 

third and fourth terms arise from the unsteady disturbances 


ahead o£ the shock which are generated by the first blade's 


motion. 





The components behind the shock are found by noting 


that the tangential component does not change 


Ve = Ve (3.4.24) 


and by using (3.4.6,19,21,22) to obtain 





Ws Stu SIN Wo (1 + = Ex 
T ES U cos%o sın?% a- = wit (314725) 
n A 9) 
+ =u- I GC Xu, SIN Yo + Y, cos ес" 


where 


Muz M SIN Ya (3.4.26) 


It is convenient to once again return to the dimensionless 
x-y velocity components | апа n , obtaining from the 


unsteady terms of (3.4.24,25) 


d 
C= my TE + SY mau + ma. V, (3.4.27) 
MU у mum Inu + NV (3.4.28) 
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where the factor e has been dropped and 
m,= — = SIN 2%, SIN? Yo (3.4.29) 
м, = er nu SIN? X, (3.4.30) 
m= cos? vo + —— ым 0 ¿E qa (3.4.31) 
ma = SER a mine a SC ЗЛО, 
M = — se (cos 2x. + Tob SIN? V, (3.4.33) 
n= RU q сот SIS, (3.4.34) 
E, (3.4.35) 
па = замма + З cost x, is 2, (3.4.36) 


Note that the above coefficients are functions of the 
local shock angle Be and so are not constant along the 
shock. 

If no unsteady disturbances exist upstream of the shock 


(u. = v. = 0) then (3.4.27,28) reduce to Teipel's single- 


I 1 
blade results. If, in addition, the shock does not oscillate 
(g = dg/dy = 0) „ then у апа v are identically zero 
indicating that there are no unsteady disturbances behind 


the shock. Equations (3.4.27,28) represent an extension of 
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Teipel's work from the single blade to the cascade, and 
they were first obtained by Chadwick for a cascade of 
wedges. It is noted for completeness that the above 
equations would apply equally well to an upper-surface 


shock wave if the signs of ту, My, N, апа п уеге 


3 
reversed. 


3.4.6 The Pressure Jump Conditions Across the Blade-Two 


Shock. An additional extension of Teipel's work was accom- 
plished by the author in order to develop a more accurate 
method of determining the flow-variable discontinuities 
than that employed by either Chadwick or Teipel. The 
improved method is treated in more detail in section 3.8, 
but a is based on the extension of the Rankine-Hugoniot 
ure relation (3.4.8) to the cascade via ‘equations 

АТ 720,21), resulting in 


^ 
P _ d2 | 
P Pay * uu -= pau a WM (3.4.37) 


where 





2.8 2. 

Te 3.4.38 

Р. d Mn SIN 2%, ( ) 
4-Ұ 2. 

ә —H— 3.4.39) 

Рг 8 + | NAUES | 

AY 2 

S (3.4.40) 
ze 
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< 2 
Pas yop Mu <от х, (3.4.41) 


Equation (3.4.37) relates the pressure jump across the 
shock to the shock displacement and the upstream unsteady 
velocities. In the single-blade case а, апа уу аге 
zero, and Teipel's result is recovered. 


Bei The Boundary Conditions for the Shock Equations. 


In section 3.8, the flow variables behind the shock are 
obtained by the simultaneous solution of the difference 
equations resulting from (3.3.14) and the shock equations 
ОООО? 28,37). In order to start the solution process, 
then, several quantities must initially be known at the 
point of intersection of the shock and ET i.e. at 


(х,у). Erom figure Dou get is.cledrethat 


I — os (3.4.42) 


Also, immediately behind the shock the flow tangency 


condition for the second blade must hold 
Y =-[1 + ce co 7 er at (Xo, Jo) (3.4.43) 


With v known equation (3.4.28) provides 


SÉ (m Ryu мам) at Xe, go) (3.4.44) 


da 
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and (3.4.44) used with (3.4.27,37) gives the remaining 


boundary conditions 


T Za" mina REC RENI mna y 


vu V "T" AR Ko, Jo) (3.4,45) 
A 
МЕ л, лс» Fits 4,4 Ра Ана у 
KS "ni SS nu | (3.4.46) 


at (Xo, ge) 


ENNcocurse for equations (3.4.43) through (3.4.46) the hatted 
quantities are just behind the shock at (колу 
Equation (3.4.45) was derived by Chadwick while (3.4.46) 


is the author's result. 


3.4.8 The Velocity Potential Jump Conditions Across the 
Blade-Two Shock.  Landahl [26] shows that the total velocity 


potential (3.3.1) is continuous across the leading-edge 
shock, but this does not imply that either the steady or 
unsteady component is individually continuous.  Teipel and 


Chadwick choose to assume that, at the first shock point 


^ 
ф.о, о) = OG go) (3.4.47) 
The author, on the other hand, chooses to find the initial 


velocity potential jump at the point (х,у) . Бог а 


perfect gas with constant specific heats, it can be shown that 
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a 


A 
P 
$ = (dir Ves. 


xn (3.4.48) 


^ ~ 

and since U, апа Р/Р. are known at (х Уд) from 
Е м 

equatıons (3.4.45,46), then so is Йй, . Figure 3.4.6 is aà 
plot of computer data for cascade T, and it shows the 
discontinuity in B4 along the shock as a function of 
distance from the blade. Note that the initial discon- 
tinuity is small justifying Teipel's assumption (325074798 
but equation (3.4.48) eliminates having to make that 


assumption at all. 


3.5 Constructing the Characteristic Grid. Because the 
characteristic directions are the same for both the steady 
and unsteady problem, the characteristic grid can be con- 
structed and employed to find both the steady and unsteady 
velocities. The grids above blade one (zone A) and above 
blade two (zone C) are composed of two families of straight, 
parallel characteristics. The mesh below blade two (zone B), 


however, is composed of one straight, divergent family ye 


and one curved family (A). 


Bese The Straight characteristic Grids Above Blades One 


and Two. Since the upper blade surfaces generate no steady 


disturbances, it is seen from (3.3.8) that, in zones A and C 


NA = хе = ME! =p" (3.5.1) 
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Figure 3.4.6 The discontinuous jump in the unsteady 
velocity potential across the leading- 


edge shock. 
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and so the characteristics in those zones are simply Mach 

waves (figure 3.5.1). Each blade has 21 equispaced upper- 
surface gridpoints, and the program requires that the leading 
edge of the second blade be positioned at one of the straight- 
mesh gridpoints. This is accomplished via the input parameters 


K and N which fix the cascade blade spacing according to 


< 
Se (32982) 
бт =ч» д + 0.03 N (3:523) 


For cascade T, K and N are 16 and 7 respectively. 
Figure 3.5.1 also shows the characteristic numbering conven- 
tion employed in the program, so ак а д ареш г may De 


located by two characteristic coordinates 


(Ly) Ta) (3.5.4) 


Note that the straight mesh is computed as if both blades 
were flat plates. The x-y coordinates of the straight 


mesh are referred to as Xl, Yl in the computer program. 


3.5.2 The Curved Characteristic Grid Below Blade Two. 
iezone B (figure 3.5.2) the characteristic grid is non- 
linear and must be constructed iteratively. Figure SV 
reveals the numbering convention employed for the curved 


mesh. The curved-mesh x-y coordinates are referred to as 
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Figure 3.5.2 The curved characteristic grid. 
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shock 
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Figure 3.5.3 Characteristic lines through the shock. 
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X2,Y2 in the program to distinguish them from those of 


the straight mesh. Before proceeding with the grid con- 
sturction, some useful relations are obtained. 


For an & characteristic passing through the shock 


from zone A to zone B (figure 3.5.3) 


A 3 
Na — Ha = Ng — Ag (375 5) 
and since A 15 zero 
3 
Nam IM? + As) (3.5.6) 
Now consider two points on the same А characteristic, 


Ecol lustrated in figure 3.5.3. From the & compatibility 


relation 


3 
NID- p Ст з= Ne (5957) 
ACL KY = ик) а А (3.5.8) 


and from the 5 compatibility relation 
3⁄2 3⁄, 
AUT, di + BIL, = ACIL K) + PAIK) (3550) 
and so it must be that 


METS AK) ALDEA (3.5.10) 


96 








which proves, in addition, that the P family is composed 
Of straight lines. 

The steady shock angle at any point along the shock is 
found by averaging the angles of the characteristic lines 


ore and after the shock, i.e. 
%= p | TAN" (X 4) e TAN" ORTA] (ИБ, М) 


Taking the tangent of (3.5.11) and applying the sum and 
double-angle trigonometric formulae leads to a quadratic 


equation in tan Ж 


R'a =o (3.5.12) 


(Ne +) TAN “2, + 2 (N iz aan, TAN X, - (A + 
The solution of (3.5.12) is rather cumbersome and is 


adquately approximated by the simpler expression 


2. 
TAN Z, = AE a A (355713) 
In the computational implementation of (3.5.13), referring 
to figure 3.5.3 again, the characteristic angles at adjacent 
points behind the shock are first averaged together, and 


the result is averaged with the zone-A angle, i.e. 


4 
2% + AMI AC, KS 


TAN ne = 


(3.5.14) 














3.5.2.1 Computing the Point (1,1). The X2-Y2 coord- 


inates of (1,1) are of course those of the blade-two 
leading edge. fe bist) is found from the steady boundary 
Monaition (3.3.5) and equation (3.3.9) while A(1,1) follows 


meen (3.5.6). 


Ec Computing the Points (1,2) and (2,2). The 


points (1,2) and (2,2) are established iteratively as 

follows. Program subroutines employed are indicated in 

parentheses, and figure 3.5.4 is provided for reference. 
1. Establish the reference point ХА on the blade, 


close to the leading edge 


ХА = Xx2 (11) + 0.005 


2. Compute A at XA , which is equal to /A(1,2) , 
using the steady boundary condition (3.3.5) and 329). 

3. Compute A at XA , which is equal to ОЕ 
using equation (3.5.6). 

4.  (Subroutine SOLVE) Locate (1,2) ‘at the intersection 
of a line through (XA, Y2(1,1)) with slope 


-\ 
^O, 2) 


and a line through  (X2(1,1), IA ES ope 


-4 
27у + АО) + A (1,222 
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то1о 
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x 


Computing curved-mesh points (1,2) and (2,2). 


(XA, Y2(1,1) ) 2082) 
a, DA С 


blade two 
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shock wave 
NO 


zone А 
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х2 


Computing subsequent shock-front points (1,1). 





Figure 3.5.4 Computing the shock-front points. 
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b. (Subroutine ITER) Guess at ЕТЕ ооб Е ТОПКО (2,2) 


em che blade 


X2 (2,2Y = X2 (VÚ 


6. (Subroutine ITER) Increment that guess 
Х2(2,2)у= %2(22) + 0001 


7: (Subroutine ITER) Compute P from 5.352559). 
8. (Subroutine ITER) Compute A(2,2) from (3.5.6). 
9. (Subroutine iTER) Solve for X2(2,2) from the 


slope БЕШЕНЕ CC characteristic through (1,2) and (2,2) 


x2(2,22-x20,2) _ 2. 
12(2,2-Y20,2 0,23% + A(2,2)% 
MO (Subroutine ITER) Does X2(2,2) = X2(2,2)? 


VES- SEO PO 11; NO: -GO TO 6. 


qd. -STOP 


ЕРО оправі па Subseguent Shock Front Points. 


Subsequent points (1,1) along the shock are also found 
iteratively. The technique follows, and figure 3.5.4 
applies. 


l. Relocate ХА on the blade 


ХА = XA t DXA 








where DXA is based on the first position of XA 


x2 (2,2) -— XA 
D%A. = ——— ы-ы 
15 


2. (Subroutine SHOKPT) Guess  A(1,I) 
AUT) = A(1, T=5 

3. (Subroutine SHOKPT) Increment that guess. 
AGL) = AUT) + 0.001 


ae (Subroutine SHOKPT) Compute EEN using ps and 
À at Pu which are the same respectively as пос 


and A(2,2) 


Аат = а MG, 2) JA 2) 


S (Subroutine SOLVE) Locate (1,1) at the intersection 


pne through (XA, Y2(1,1)) with slope 


mj 
A(1, 1) 


ЧОЛО е through  (X2(1,1-1),Y2(1,i-1)) with slope 


-4 


ZN AO I-OU + AU, Ty 
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6. (Subroutine SHOKPT) Compute A1, D From che 
steady boundary condition (3.3.5) at XA and equation 
E 35.9). 

7. (Subroutine SHOKPT) Does fat) . 
m5: GO TO 8; NO: GO TO 3. 


ОО STOP. 


55.2.4 Computing the Curved-Mesh Interior Points. A 


point (I,J) interior to the curved mesh is located (figure 
E»)-at the intersection of a line passing through 


(I-1,J) with slope 


| 
A (I-1,3)4 


БО а line through (1,1) with slope 


| 
A E) 


Po ean Computing Subsequent Blade-Surtace Points, 


Subsequent points of the form (I,I) on the lower surface 
are computed in a manner exactly analagous to the computation 
of the point (2,2) as described in steps 5 through 11 of 


section 3.5.2.2. 


КО ОШ О ШЕ по the Steady Velocities and Pressures. With 
/^ and A known throughout the grid, the steady velocities 


BNET lir Ci Ty from (3.3.8,9). The steady pressures on 


(Res 








(1-1,1-1) (1,1) 


blade two 


(I-1,J) 


B 





higusess.5.5 Computing the curved-mesh 
ето еро 


the upper surface of blade two are zero and on the lower 


surface are given at the surface points  (I,I) by 


Cp = - 2 ио ст) (3. 631) 


° 


Computing the Unsteady Flow Variables at the 
Straight-Hesh Grid Points. 


3.7.1 The Blade-One Upper-Surface Grid Points. In the 


straight mesh the unsteady velocities and velocity potential 


are designated by 


u = UllR 4 i UllI Vi = VllR + 1 V11I 


б =z FÌR 4 i FlI (3. 7251) 





The unsteady boundary condition for blade CIO A 


allows determination of the MI velocities 


VI1R(1,1) = -1  VLLI(I,I) = -k|x1(1,1) - b] (36772) 


Using the А compatibility relation (3.3.14) in difference 
form results in a system of linear equations for the 

EnEnowns UllR(I,I), UllI(I,I) . Figure 3.7.1 A applies, 
and the equations appear in matrix form in section A.l of 

the appendix. The velocity potential is found by integration 
of the velocities, using equation (3.3.18); see section A.2 


of appendix A. 


3.7.2 The Straight-Mesh Field Points. As seen in figure 


DOS, grid points not on the blade require use of both 
compatibility relations. The equations resulting for the 


unknown u and v velocities at each field point (I,J) 


1 l 
are given in A.3 of the appendix. Equations (3.3.17,18) 
are averaged to find the velocity potential; see section 


A.4 of the appendix. 


3.7.3 The Grid Points Immediately Upstream of the Shock. 
Figure 3.7.1 C shows that the shock points lie in the straight 
mesh and therefore the flow variables at each shock point — 
upstream of the shock — are found by means of an interpolation 
algorithm. The distances RI in figure 3.7.1 C are computed 
and any desired quantity Q(1,K) is found from the values 


EOS at the four surrounding points 


NA 






U3 ,V4 unknown 


Ulv} known 


V 


(I-1,T 
P uj unknown 


v. known from 
boundary cond. 


blade one 





(2,1) 





blade two 






Blade two E 


ex 
EXCI-l,J) 










flow veriables 

here are about 

equal to those 
at leading 


edge behind 
N shock, 
e 
here are about 
equal to 


A 
those at e VI) 
(L,I) behind 
the shock. (1,141) 


shock 


flow 
variables 





D 


Figure 3.7.1 Computing the straight-mesh and 
curved-mesh grid points. 
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RRLQUI-13-D+RR2 Q(T-,T)+RR3 Q(I,J)+RRAY QUI,I-1) 


A ee Ee 
Ber E (3.7.3) 


where 


4- 
RRS-TT RI Ro-T.RRJO (3.9.4) 
I+J Sl 


3.8 Determining the Jumps in the Flow Variables Across the 
Blade-Two Shock. It is in the determination of the jump 
conditions that the author employs a technique — alluded to 
in section 3.4.8 — that is different from the method used 
Eo Teipel in the single-blade case. The essence of Teipel's 
approach is to: | | 

l. Assume that the velocity potential does not change 
Beross the shock at the first point (3.4.47). 

2. Use the shock velocity equations (3.4.27,28), with 
boundary conditions (3.4.42,43,45), in conjunction with the 
A compatibility relation (3.3.14), to solve for Ж SR 
(behind the shock) and g at each shock point. 

3. Use equation (3.3.18) behind the shock to recover 
ГД .at each shock point. 

In contrast, the method employed in the present study 1s tos 


A : 
1. Use (3.4.46,48) to determine Йй, behind the shock 


de the first shock point. 
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2. Use the shock velocity equations (3.8 221,28)... the 
shock pressure equation (3.4.37), the P compatibility 
relation (3.3.14) and the boundary conditions (3 4.42. 43; 
45,46) to solve for ЕУ g and P/P at each shock point. 

3. Use equation (3.4.48) to find A from the pressure 
jump at each shock point. 

The latter method results in a larger system of equations 
at each point, but no assumptions concerning the continuity 


of the velocity potential are required. 


ENS Details of the Jump-Condition Computations. 


Step 2 above requires some additional comment. The shock 
equations (3.4.27,28,37) are applied in difference form from 
one shock point to the next. They are the last six equations 
in the matrix equation given in appendix section A.5. The 
R relation supplies the first two equations in that matrix 
equation. Note from figure 3.7.1 D that using the A relation 
at shock point (1,1) requires two tacit assumptions. First, 
the flow variables at the point where the A characteristic 
through (1,1) leaves the blade are essentially those at 
the point (1,1) immediately behind the shock. Second, the 
flow variables on the A characteristic through (1,1+1) 
and just behind (1,I) are about the same as those at 
(1,I) behind the shock. 

In the computer program the flow variables at shock 


point (1,I) , ahead of the shock, have the form 


етот т) ота т) (3.8.1) 





while, at the same point immediately behind the shock, 


they are designated 


ur = URJCI) + i UIJ(I) Cece) 


MEL stored in the form of (3.8.1). 


3.9 Computing the Unsteady Flow Variables at the Curved- 


Mesh Grid Points. 

3.9.1 The Blade-Two Lower-Surface Grid Points. The 
unsteady boundary condition for the second blade provides 
the lower-surface Vi velocities (figure 3.9.1 A). The 
equations are given in section A.6 of the appendix. The 


u velocities are found via the & compatibility relation 


el 
(section A.7) and the velocity potential follows from 


eam@ation (3.3.17) (appendix section A.8). 


3.9.2 The Curved-Mesh Field Points. Both compatibility 
relations are used at a field point, as seen from figure 
3.9.1 B. The resultant matrix equation appears in appendix 
section A.9. The velocity potential is found from equations 
(3.3.17,18) (appendix section A.10). 

With the curved mesh completely computed the unsteady 
pressures at blade-two, lower-surface grid points are given 


by 


Ee -2 [UL2R(1,1) - k F21C1,1)] 


(3.9.1) 
25 |UL21(1,1) + k F2R(1,1)]| 
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Figure 3.9.1 Computing the curved-mesh grid points. 


3.10 Results and Conclusions. 

3.10.1 The Isolated-Blade Case. If the unsteady boundary 
condition along the first blade is set to zero, the XO 
blade sees freestream conditions and its unsteady pressures 
can be compared with Teipel's isolated-blade results [23]. 
This is accomplished in figure 3.10.1, where it is noted 
that a small discrepancy between the two solutions exists. 
This discrepancy is attributable in part to the different 
techniques used by Teipel and the author in solving the 
shock equations (see section 3.8). But it is also noted 
that Chadwick [5] too encountered a discrepancy when he 


applied Teipel's method to the single,. biconvex airfoil. 
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Figure 51071: Isolated-Blade Results. 


Teipel points out that; for Пе cue plade, as К 
gets ge the linear and nonlinear solutions approach 
each other. This phenomenon is also evident for the second 
blade in cascade T, a5 сееп 1П figure 331022 which compares 


Platzer's linear solution [17] with the nonlinear results of 


this inyestigation fora K of 1.5 - 


3.10.2 The Stead Pressure Results. Figure 3.10.3 shows 


the steady pressures for the lower surface of blade two. 


It is seen he agreement between the experimental results 








Cascade T te 1.55 


ie 


blade 2 lower surface 


ES 1.5 J: 30.09 


Platzer-P Strada-S 





Figure 3.10.2 High-k comparison of the linear 
and non-linear solutions. 


Of Fleeter [25] and the numerical results of the computer 
model is -good Beer RU fair at the leading edge. 
The program results compare favorably with theoretical 
predictions from second-order Prandtl-Busemann theory and 
a third-order theory due to Ferri [30] which takes entropy 
losses into account. This latter comparison again demon- 
Strates that the computer solution is not limited by its 
use of the isentropic Rankine-Hugoniot pressure relation 


(4.8). 


3.10.3 The Unsteady Pressure Results. Figures 3.10.4 
ООО j j 3 10.12 show the computer results of the present 


study, the experimental data from Fleeter's investigation, 
and the linear results of Platzer's method-of-characteristics 


program for nine interblade phase angles. The pressure 











MO T: 0,03 
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theory Ea 


Fleeter © Strada (a 





Figure 3.10.3 Blade-two lower-surface 
steady pressures. 


data are presented as an amplitude and phase angle relative 
to the motion of the second blade. The phase angle is 
measured from the point where the second blade attains 
maximum up pitch to the point of peak blade pressure. Along 
the upper surface of blade two the linear and nonlinear 


numerical results agree exactly — because the upper surface 


TIZ 








is flat — providing an additional check of the nonlinear 
program. Along the lower surface the nonlinear effects 

show up as differences between the two predicted pressure 
distributions — differences which vary with interblade phase 
angle. The predicted nonlinear pressure amplitude is, in 
all cases examined, at least as large as the linear predic- 
tions, which is as expected. In figure 3.10.5 (A= -4.4°) 
the linear and nonlinear amplitude predictions are nearly 
coincident, while in figure 3.10.9 (A= 130.6°) they differ 
considerably. Figure 3.10.13 is a plot of the difference 
between the nonlinear and linear amplitude predictions versus 
the turning angles experienced by the flow at the leading 
edge of blade two. It shows that the nonlinear computer 
predictions .dre consistent since the variation from the 
linear solution ос es with the flow turning angle. 


That turning angle is given by 
TA.= 57.3% (1- COSA) + 6.8" (3.10%) 


when blade two is in the maximum up-pitch position. Note 
that the large turning angles in figure 3.10.13 result from 
the choice of a one radian maximum pitch amplitude for 
nondimensionalization purposes. 

Comparison of the nonlinear computer predictions with 
Fleeter's experimental data, however, reveals that the 


theory and the experimental results do not consistently agree. 
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Figure 3.10.4 Blade-two unsteady pressures. 
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Figure 3.10.5 Blade-two unsteady pressures. 
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Figure 3.10.6 Blade-two unsteady pressures. 
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Figure 3.10.7  Blade-two unsteady pressures. 
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Blade-two unsteady pressures. 
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Figure 3.10.9 Blade-two unsteady pressures. 
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Figure 3.10.10 Blade-two unsteady pressures. 
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Figure 3.10.11 Blade-two unsteady pressures. 
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Figure 3.10.12 Blade-two unsteady pressures. 
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Figure 3.10.13 Nonlinear variations from linear theory. 


In figure 3.10.11, for example, there is adequate upper- 
surface agreement with the theory, but the lower-surface 
pressure-amplitude comparison is poor. In fact, in no case 
is there agreement between data and theory for both phase 
and amplitude along the lower surface. The lower-surface 
experimental data appear to be somewhat inconsistent. They 
sometimes vary widely from one point to the next for the 


same interblade phase angle (figure 3.10.6) and they also 
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vary a great deal for rather small changes in interblade 
phase angle (figures 3.10.10,11,12). In addition, in each 

of the figures there is at least one amplitude data point 

for the nonlinear (lower) surface that is below the corres- 
ponding linear-surface data point. This phenomenon is 
contrary to expectation because the thickness surface must 
generate a larger magnitude disturbance than the flat surface. 
Certainly, additional experimental data are required to 


adequately evaluate the theoretical methods. 


3.10.4 Conclusions. In summary, it is concluded that the 
extension of Teipel's single-blade methods to the nonlinear 
cascade results in: 

l. Steady pressure predictions which are in accord 
with К лег nonlinear theories and which are supported, to 
some extent, by experimental results. 

2. Unsteady pressure predictions which exhibit expected 
trends relative to linear theory, but which are not consis- 
tently in agreement with experimental data presently available. 
The significant influence of thickness on the unsteady 
pressures is apparent in both the theory and data. 

3. A numerical approach to the unsteady cascade which 
appears to have the potential to treat the effects of blade 
thickness, blade camber, and steady angle-of-attack loading 
to provide a more realistic model of the fan and compressor 


blades in a turbomachine. 
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APPENDIX A 


A.l Equations for the Unknown uj Velocity at Blade-One 


Surface Points. 
E NPererences: sectdon 5.7.1: ficure 3.7.1028 


2, AUxil®ary constants: 


AL = 2kM2(0.025)/A,  A2 z k“M“(0.025)/Aa 
\ 

A5 z 190.0125 A2 AG = АЗ/М2 

А5 = 1— 0.0125 A2 A6 = AS5/AR 


3. Matrix equation: 
A ira, ERT EE 
A1/2 AS ULLI(I,I)| {A3 ULLICI-1,1I) = 
Ai ULIICI-L,1)/2 = A4 VLIR(I,1) +A6 VLIR(I-1,1)+ 


mead) Whee (I-1,1)/2 = A4 VILICI,1) *A6 VLLECI-1,1)+ 


+Ф А2 Р1Р(1-1,1) 


42 FlI(I-1,1) 


A.2 Equations for f4 at Blade-One Surface Points. 
1. Reference: section 3.7.1; figure 3.7.1 A. 
2. Equations: 


FlR(I,I) 


FIRCI-1,1) + {ULIRCI,T) + ULIR(I-1,1) - 
- [vl1R(I,I) + vl1R(I-1,1)] AR (©0125; 
miei, =) S PLICTSL, I) + {ullıcı,7) 4 ULLICI-1,1) - 
EDT CEU E anro | Pre (0.0125) 
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A.5 Equations for the Velocities uy and vj; at Straight-Mesh 


Grid Points. 
1. References: section 3.7.2; figure 3.7.1 B. 


2. Auxiliary constants: 


Al z 2kM^ (0.025)/A А2 = k2M2(0.025)/A 

A5 = 1 — 0.0125 A2 А& = 1 + 0.0125 A2 

AS = А5/ ^2 Аб = Ae fla 

5. Matrix equation: 

A5 J -Al/2 -Аб 0 ULIR(I,J) AG UllR(I,J-1) + 
AL/2 АЗ О -A6 UL1I(1,J)|  [A4 UllI(I,J-1) - 
А5_ -А1/2_ 46 O EE RT yee 
А1 /2 АЗ О Аб VUE TD A4 ULLICI-1,J) ~ 


«* AL ULlI(I,J-1) — A5 VIlR(I,J-l1) + A2 FIR(I,J-1) 
Ф 

- Al ULl1R(I,J-1) — A5 VllI(I,J-1) - A2 FlI(I,J-1) 
2 

+ AL ULLICI-1,J) + A5 VLIR(I-l,J) + A2 FlR(I-1,J) 
2. 


~- Al ULIR(I-1,J) + AS VILICI-1,J) = A2 FLICI-1,J) 
2 


eea Equations for gu at Straight-Mesh Grid Points. 
FIR(I,J) = $[FIR(I,J-1) + FIR(I-1,J)| + 0.0125 ULIR(I,J) + 
+ 0.00625 [ULLR(I,J-L) + ULIR(I-1,3)| + 
+ 0.00625 [V1IR(I,J-1) - VLIRC(1-1,J)] fx 
FLI(I,J) = 2 [F11(1,3-1) + FLI(I-1,3)] + 0.0125 ULLICI,I) + 
о оосар аата e 
4 0.00625 [111(1,7-1) - v111(1-1,3)]/9% 


1. Peferences: section 3.7.2; figure 3.7.1 B. 
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A.5 Equations for the Jump Conditions across the Blade-Two 
Leading-Edge Shock. 
1. References: section 3.8.1; figure 3.7.1D. 


2. Auxiliary constants: 


Y Y 
2/ [Ma 1) + A1,1-1)) 


AL = 

A2 2 4? [x2a, - x20,1231)]/ [a, D « ^a, 1-1)] 
BENZ д2 к/2 

py = [y2c1,1) - y2c1,1-1)]/2 


3. Matrix equations: 


1 mr No А572 Q О 0 URJ(I) 
A2/2 1 O Al O -A3/2 0 О UIJ(I) 

1 0 OO 0 -m/DY m,  - ||VRJ(I) 

0 1 dO 6 o ma -mj/DY||vixn|. 
0 0 NEC 0 0 -nj/DY п, PRICI) | 
О 0 O 0 -n5 -ni/DY| | F1J(1) 
ҳм2 о 0 0 0 $M^k -ру/рү р, GRCI) 

0 -4M2 о о Mk 0 -P, -p, /DY| | GI(I) 


URJ(I-1) + A2 UIJ(I-1)/2 « Al VRJ(I-l) + АЗ ЕВЈ(1-1)/2 . 










UIJ(I-1) - A2 URJ(I-1)/2 a Al VIJ(I-1) + A3 FIJ(I-1)/2 

m3 [UL2R (1,1) + UL2R(1,1-1)] + m¿[v12R(1,1) + V12R(1,1-1))- 
m3 [UL21(1,1) + uL21(1,1-1)] + ma [у121(1,1) + v121(1,1-1)] - 
ma [UL2R (1,1) + UL2R(1,1-1)] « n4 [v1220, D + VI2R(1, I-1)] - 
my [UL21(1,1) + 0121(1,1-1)] + ny [v121(1,D) + v121(1,1-1)]- 
(р, - м2)[112в(1‚1) + UL2R(1,1-1)] + py[v12R(1,1) + 


(ps - M2)[U12IG1,I) + 0121(1,1-1)]) + ра [У12101,1) + 
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A22 (Continued): 
3. Matrix equation (right hand side continued): 





URJ(I-1) 


ту GR(I-L)/DY - m, GI(I-1) 


UIJ(I-l) - mj GI(CI-1)/DY 4 m4 GR(I-1) 


VRJ(OI-1) - n, GR(I-1)/DY - nj GI(I-l) 


VIJ(OI-1) - ni GI(I-1)/DY - no GR(I-1) 
VL2R(1-1)] + kYm? [F21(1) + F21(1-1)] +4 MURICI-1) = 


> 
+ V121(I-1)| - k$M? [ron c1) + F2R(1-1 )]* x M^UIJ(I-1) 4 








"Ж 





-$M^k FIJ(OI-1) - ру GRCI-L)/DY - p, GICI-1) 


+%#M°?k FRJ(I-1) - p) GI(I-1)/DY + р„ GR(I-L) 





A.6 Boundary-Condition Equations for the vj Velocity at the 
Blade-Two Surface Points. 
l. References: section 3.9.1; figure 3.9.1 A. 
2. Auxiliary constants: es. CSC EI 
3. Equations: 
VL2R(1,1) = - COS(M) - k(X-D)SIN(A) 
v121(1,1) = - SINGA) = k (X-b)C0S (41) 
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A.7 Equations for the uj Velocity at Blade-Two Surface 


Points, 

L. References: section 3.9.1; figure 3.9.1 A. 

2. Auxiliary constants (also used in A.9 equations): 
DXA = X2(1,J) - X2(1-1,J)  XLa= 3[M1,3) + A1-1,3) 
DXB = Х2(1,Ј) - X2(1,J-L) XIB = 3[M1,J) + M1,3-1)] 
VAR = 3[y12R(1,J) + VL2R(1-1,3)] пв = №І,Ј) - AG,3-1) 
vat = 2[у121(1,7) + VL21(1,3-1)] AL = 2 км2 
DVAR = VL2R(I1,J) - VL2R(I-1,J) A2 2 k Al/2 
ОЧАТ = VL2I(1,J) - V121(I-1,J) 
SLA = 4 [ca Aa- A] ans Afen af ausa] 
р1хв = [А‹(т,л+1) - мт, о] Бост, ож) Eco 


3. Matrix equation: 


2 | S 
DLXB DXA . A2 DXA _ Al DXA 
re Эхх LA UL2R(1,J) 
Al DXA DLXB DXA A2 DXA? 
ШО CTA l + eelere SC UL2I(1,J) 





LXB DYA A2 XA? 
UL2R(I jp A 5 5 1+ 12161-1, Ј) = as 


2 ` 
aart pit - DLXB DXA 4 A2 DXA |- 1 _1. ту А1 ЮХА 








DVAR _ A2 VAR DXA2 A2 DXA 
+ SIR ta 4 Р2В(І-1,Ј) AG 


2 
DVAI A2 VAI DXA A2 DXA 
УСТА Y XLA SLA + F21(1-1,J) VLA 





4. Note: at surface points I = J. 


129 








A.8 Equations for @)_at Blade-Two Surface Points. 
1. References: section 3.9.1; figure 3.9.1 A. 
2. Auxiliary constants: as defined in A.7. 
3. Equations: 
F2R(1,3) = F2R(1-1,3) + 44[UL2R(1,3) + UL2RCT-L, 3)] 
+ VAR/SLAY DXA 
BAR, I) = F2I(I-1,J) + (3 [u21(1, 9 + UL21(1-1,J)) 
+ VAI/SLAY DXA 


4. Note: at surface points I = J. 


A.9 Equations for the Velocities uy and V4 at Curved-Mesh 


Grid Points, 


1. References: section 3.9.2; figure 3.9.1 B. 
2. Auxillary constants: as defined in A.7; 


3, Matrix equation: 






l 4-.DLXB DXA ... A2 DXA? = LAL DXA 1 __А2 DXA? 
4 XLA 2 XLA 2 XLA SLA 2 XLA SLA 


Al ЮХА у ¿ DLXB_DXA _ А2 DXA? 0 Y __A2 DXA? 

2 XLA 4 XLA DEREN SLA 2 XLA SLA 

E A EE 0 
ZEE 2 XLB 2 XLB SLB 2 XLB SLB 

Al DXB p DB A2 DXB? А 14. А2 рхв2 

2 XLB 4 XLE 2 XLE SLB ' 2 XLB SLB 


(Note: continued on next page.) 
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A.9 (Continued): 


3. Matrix equation (unknown vector and right-hand-side 


matrix continued): 


NL2R(I,J) Ul2R(I-1,J) h .-PLXB DXA 1 к > 


4 XLA 2 XLA 














oan, |uarcr-1,3) [1 - 25 XA y 42 002 
Vl2R(CI,J) 0128(1,Ј-1) Ë Ë TUE «ALES k 
eT (ie) UL21(1,J-1) [1 Я "ux « A38 | 

+ U121(1-1,J) LDL + VL2R(1-1,3) к 
user 5) AL DZA y VIKC- 1,3) - = SLA + т - 


~~ UL21(1, 5-1) AL DXB 4. vl2R(I, J-1) [+ — A2 DXB2 | 4 
2 XLB SEB 2 XLB SLB 


< 


2 
— UL2R -1) AL DXB + v121(1,J-1) —_A2 DXB + 
B---5 E : TI — 2 XLB SLB 














4 F2R(I-1,J) AZ DA 
XLA 


ZA 


+ F21(1-1,3) A2. DXÀ 


4% 


4 F2R(I,J-1) .à2 DXB 
XLB 


F21(1,J-1) A2. DXP 
Е : XLB 
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A.lO Equations for Dy at Curved-Mesh Grid Points, 


і. References: section 5.9.2; figure 5.9.1 B. 
2. Auxiliary constants: as defined in A.7. 
5. Equations? 
F2R(1,3) = 3[F2R(1-1,3) + F2R(1,3-1)] 
hs San zc, EE I DO + YLZRCL,I) + у12в(1-1, ee, 


+ HUL2R(1,J) + ЇЛ2Р(Т,1°10) = о e 


F21(1,3) = 3[F21(1-1,3) + F21(1,3-1)] 


К FUL21(1,3) + UL21(1-1,J) + YL21(1,D) 4c VI2ICI-L J| pxa 


^ 4{U121(1,9) + UL2I(I,J-1) = EIC aa uerb ox 
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